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Abstract 

Using a three-level algorithm we perform a high-precision lattice computation of the 
static force up to Ifm in the 2+1 dimensional SU(5) gauge theory. Discretization errors 
and the continuum limit are discussed in detail. By comparison with existing SU(2) and 
SU(3) data it is found that otq = 1.65 — ^ holds at an accuracy of 1% for all > 2, 
where tq is the Sommer reference scale. The effective central charge c(r) is obtained 
and an intermediate distance rg is defined via the property c(rs) = It separates 
in a natural way the short-distance regime governed by perturbation theory from the 
long-distance regime described by an effective string theory. The ratio rg/ro decreases 
significantly from SU(2) to SU(3) to SU(5), where rg < tq. We give a preliminary 
estimate of its value in the large- limit. The static force in the smallest representation 
of AA-ality 2, which tends to the k = 2 string tension as r — > oo, is also computed up 
to 0.7fm. The deviation from Casimir scaling is positive and grows from 0.1% to 1% in 
that range. 



1 Introduction 



One of the prime observables giving an insight into the workings of four-dimensional 
SU(A^) gauge theories is the force F{r) which two static quarks separated by a distance 
r exert on each other. It may for instance serve to define a renormahzed couphng 
constant, and hence the knowledge of F(r) at all distances allows one to extract in 
units of the string tension a, establishing the connection between the regime accurately 
described by perturbation theory and the regime thought to be described by an effective 
string theory. While in practice finite-volume schemes are better suited to overcome the 
hierarchy problem in connecting the very-short-distance regime to the non-perturbative 
regime F(r) is ideally suited to probe the string properties of the chromo-electric 
flux lines, since its functional form at asymptotic distances is a central prediction of 
effective string theories O IHj . 

Although the physically most important SU(A^) gauge theories are undoubtedly de- 
fined in d = 4 dimensions, in this paper we shall focus on the d = 3 theories. Apart from 
being interesting in their own right (they exhibit a mass gap and linear confinement, as 
demonstrated numerically in ^), the N = 3 case describes the magnetic sector of QCD 
at asymptotically high temperatures and is thus relevant to the description of real-world 
physics (see ^ for a review on this aspect). 

The theme of this paper is the A^-dependence of static forces. Due to super- 
renormalizability in d = 3, the quantity analogous to ^/a/Ay^ is y/a/{g'^N), which 
has been calculated in Its A^-dependence is well described by a constant, plus small 
0(1/A^^) corrections. This means that once two SU(A^) theories have been matched 
in the ultraviolet (as originally proposed by t'Hooft 0), their far infrared behaviour 
will agree too; and vice versa, these statements holding up to Oil/N"^) corrections. 
The low- lying spectra of these theories, compared in units of ^/a, are also very simi- 
lar I^IHIHl- The static force F{r) gives us an independent way to compare the theories 
at an adjustable distance scale. 

The asymptotic approach to a constant force is of particular interest, because the 
Liischer term [5], as well as the next term in the (or^)"^ expansion are universal. 
Thus by comparing the static forces of different SU(A^) gauge theories, one obtains a 
test of universality. By the same token, the pre- asymptotic differences between these 
functions of r inform us on the A^-dependence of higher order string corrections, as well 
as on the terms that vanish at long distance faster than any power of 1/r. Indeed it 
is widely believed |12| I13j that the flux-tube also admits massive modes, in particular 
longitudinal compression modes. How their effects are to be disentangled when only a 
finite number of terms of the asymptotic series are known is however not clear to us. 

What can we expect about the A^-dependence of the higher-order coefficients in 
the 1/r expansion of the static force? Since the latter is related by open-closed string 
duality [21 to the spectrum of torelons in the theory defined on a spatial Lx oo cylinder, 
and that the spectrum is expected to have a finite large- A^ limit with 1/A^^ corrections at 
all L, one would expect (although it is not mathematically guaranteed) the coefficients 
of the 1/r series to have a finite large- A^ limit with 1/A^^ corrections. Note that the 
perturbative series in g'^N, which is also asymptotic, does have this property, at least in 
the low orders where the coefficients are computed explicitly. 

Powerful Monte-Carlo techniques were developed in [14^ that triggered a very ac- 
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curate determination of F{r) in the range 0.2fm < r < 1.2fm (throughout this paper 
the 'fermi' is identified with 2ro, where rg is the Sommer reference scale jl5j). Here we 
extend these techniques to address the following physics issues: 

• what is the A^-dependence of F{r)/a? A special case is the separation r = ro, 
where this amounts to studying the A^-dependence of TqCt; 

• since the A^-dependence of TqCt turns out to be very weak, we then ask about 
the A^-dependence of the effective central charge c(r) = —^f^-^r^- Any effective 
theory whose light degrees of freedom correspond to the transverse fluctuations of 
the flux-tube predicts lim,.^oo c(r) = ^(d — 2), independently of N ^0]. However 
the approach to this asymptotic value does have an AT-dependence which turns 
out to be substantial; 

• what is the dependence of the static force on the color representation of the 
charges? Many confinement models predict the static force in a representation R to 
be proportional to its quadratic Casimir Cr ^1 E[ El El EOI ■ This question has 
been addressed extensively (in d = 4 [HI HI HSl , but also in d = 3 PH 1^ 1^ l?7j ^ : 
our aim will be to quantify the small deviations from the 'Casimir scaling' predic- 
tion. 

Concerning this last point, we focus on the representations of two 'quarks' symme- 
trized or antisymmetrized in color in the d = 3 SU(5) theory. These are particularly 
interesting, as the corresponding static forces are known |2H] to be linearly confining 
with a string tension larger than in the fundamental representation. Since perturbation 
theory predicts Casimir scaling of the static forces to hold at short distances and 
that the string tensions have been found to be near-proportional to the smallest Casimir 
of the given AA-ality [281 l3Uj . it is a priori plausible that Casimir scaling of the static 
force provides a good approximation at all separations r in the antisymmetric case. The 
situation is thus different from studies in SU(2) or SU(3), where the static force in any 
representation asymptotically vanishes or is given by the fundamental string tension. 
The difference is due to the Z{N) symmetry, which protects the k = 2 string from 
screening once N > 4. 

Note that rgcr, c(r) and F^{r)/Fp{r) (the index refers to the representation, F being 
the fundamental one) are all quantities which have a continuum limit - unlike for instance 
roV{r), the static potential in units ofrQ^. Because we are after small effects (corrections 
of c(r) from vr/24, deviation oi F^{r)/Fp{r) from Cr/Cf), the O(a^) discretization errors 
need to be at least estimated, or even better the continuum extrapolation must be carried 
out. 

Accurate SU(3) data was obtained a few years ago with a multi-level algorithm in 
d = 3 and d = 4 jlUj (in the latter case also more recently in ^T|), and SU(2) data 
is also available in d = 3 |3H I33j . Here we perform a. d = 3 SU(5) calculation that 
supplements the existing data. We do so by generalizing the multi-level algorithm [TJj 
for the higher representations (see also 1^ for the adjoint case), and introduce a further 
level of factorization |34| I35j of the Polyakov loop correlator from which the static force 
is extracted. 

The algorithmic and technical details are given in section |21 The new static force 
data for SU(5) is presented in section 13.11 and the A'"-dependence of F{r) discussed 
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in 13.21 Sections 13.31 and 13.41 concern the effective central charge; in the latter the SU(5) 
data is compared to the existing SU(2) and SU(3) data. Section |1] discusses the ratios 
of static forces and their deviations from Casimir scaling. Our conclusions are gathered 
in section |S1 The appendix contains a discussion of discretization errors and of the 
continuum limit. 



2 Lattice simulations 

In this section we describe the technical details of the computation. The number of 
dimensions d and the number of colors N are still kept unspecified at this stage. 

2.1 Action and observables 

We use the standard plaquette action: 

S = |ER«Trp{l-C/,}, /3 = ^. (1) 

plaq ^° 

where Up G SU (N) is the ordered product of links around a plaquette. The quantity 
TtrC/ is the trace of U in the representation R; the subscript R=F corresponds to 
the fundamental representation. The size of the lattice is L'^~^T, and the boundary 
conditions for the link variables are periodic in all directions. 
Our primary observable is the Polyakov loop correlator 

(^R(^)^R(y)) = 1/ D[U] P^(x)PR(y)e-^[^], D[U] = YldU,{x) (2) 

where dU is the normalized invariant measure, Z is such that (1) = 1 and 

Pk{x) = Tvr{Up{x)}, Up{x) = Uo{x)Uo{x + a6)...Uo{x + {T-a)6). (3) 

The trace may be taken in a general representation R of SU(A^). Since the link variables 
belong to the fundamental representation, we relate TrRf/ to the trace in the funda- 
mental representation. Here are some examples for three irreducible representations: 

Tta{U} = \Ttf{U}\^-1 (4) 
Tr2a{f/} = i((TVF{f/})'-TVF{f/'}) (5) 
Tt2s{U} = i((T¥F{f/})' + TYF{f/'}). (6) 



The representation A is the adjoint representation, 2a and 2s are the representations of 
two quarks (anti)symmetrized in color. We have computed the Polyakov loop correlator 
for R=F, 2a and 2s. Throughout this paper the index '1', or an index omitted altogether, 
stands for the fundamental representation F. 
We now describe the algorithmic details. 
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2.2 The multi-level algorithm 

Rather than alternating fuh-lattice sweeps with measurements, the algorithm we employ 
proceeds in a multi-level scheme which exploits the manifest locality of the action 
to factorize the Polyakov loops correlators into several more local functionals of the 
gauge field. The short distance fluctuations of the fields may then be averaged out 
independently on each of these factors. 

2.2.1 Correlator in the fundamental representation 

The factorization which yields the largest gain is the one introduced by Liischer and 
Weisz jl4(. In their algorithm, the two (untraced) Polyakov loops are sliced up into 
no = T/Aq segments by no equidistant time-slices. Recall the direct product of two 
matrices, 

{U V)a/3y5 = Uai3V-y5. (7) 

We first introduce the one-line operator 

Tl(x) = T Uf:t' Uo{x[,] + ajb). (8) 

It depends only on the temporal links between time-slice sAq/o and (s + l)Ao/a, time- 
slice T/a being identified with time-slice 0. For a point x with xq = 0, we defined 
X[s] = X + sAqO, s = 0, . . . ,no — 1; X and y can of course be taken to lie in the xq = 
time-slice without loss of generality. The symbol T means that the product is time- 
ordered. Next we introduce a two-line operator: 

Ti\x,y) = Tl(x)*®Tl(y). (9) 

Recall the product of two such direct-product matrices 

{U^V)-{U'0V') = {UU')(g){VV'), (10) 

and the trace 

Tr{U CSV} = Tt{U} Ti{V}. (11) 
The Polyakov loop operator in the fundamental representation can then be written as 

PHxyP^iy) = Tr {t UTio' Ti\x,y)} . (12) 

Because the contribution to the action of a particular link variable only depends on 
the local staples, the Boltzmann weight e""^'*^] also factorizes and one can write 

(PF(x)*PF(y)> = (Tr{T UT=o' {T-^ix,y))s}) (13) 

where {.)s denotes the average with the same Boltzmann weight but with fixed spatial 
links in time-slices sAq/u and (s -|- l)Ao/a. This formula is the basis of the algorithm 
introduced by Liischer and Weisz (14l| . 

We shall exploit an additional factorization level [341 135j . by which (T|T(x,y))s can 
be estimated as follows. Single out the direction 1 and decompose the slab of volume 
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Figure 1: The geometry of slabs and blocks in the d = 3 case. The thick lines symbolize 
the Polyakov loops, and A_Lnj_ = L. 

Ijd.-^ . Aq into ni = L/Ai smaller blocks of volume Ai • L'^~'^ ■ Aq. Fig. nmay help to 
visualize the situation. We then have 

(Tf = (Tl(x)):,^ ® {^l{y))sby. (14) 

The blocks labeled by the integers hx = [xi/Ai] and by = [yi/Ai] must be different, 
with X belonging to the former and y belonging to the latter. The expectation value 
{.)sb is taken with fixed spatial links in time slices sAq and (s + l)Ao and also fixed links 
in the set 

{U^{z) I zi = 5, • Ai and /i ^ 1} . (15) 
2.2.2 Higher representations 

To incorporate the Polyakov loop correlator in higher representations in the multi-level 
scheme described above, higher direct products must be computed. We generalize Eq. [3 
so that, for any two color tensors U and V of rank n and m, (the color indices taking 
values 1 to A^), 

{U ® V")ai,...,a„;/3i,...,/3„ = Vai,...,/3™ • (16) 

With the definition 

Tr(x,y) = T|(x)®Tl(y), (17) 
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we can define the AA-ality 2 version of tlie two-line operator 

TpT(^^y) =TlT(x,x)*®TP(y,y), (18) 
which has x entries. The product of two such objects is defined by the property 
(([/ U') {U" U'")) ■ {{V V) {V" ^ V'")) = {UV ^ U'V) ® {U"V" ^ U"'V"'). 
and multihnearity. 

The observables needed to compute the correlator in the 2a and 2s representations 
may now be expressed as 

P^xTPky) = Tr2i{Tn:£o'TPn^,y)}, (19) 
TT{UpixfrTr{Upiyf} = Tr,,, {t UT=o' T^i^^Hx,y)} , (20) 
P|(x)*TV{C/p(y)2} = Tr2iii{Tn:io'TF(^>y)}> (21) 
where the trace operations are defined such that 

Tr2i{{U0U')(S){V(S)V')} = Tv{U} Tt{U'} Tt{V} Tt{V'}, (22) 
Tr2ii{{U0U')(g){V^V')} = Tt{UU'}Ti{VV'}, (23) 
Tr2iii{(C/ U') (g){V(g) V')} = Tr {U} Tr {U'} Tt {VV}. (24) 

For the adjoint representation, one needs the definition 

TP(:E,y) =Tf (x,x)*®Tf (y,y). (25) 

Then 

mx)PF{y)\' = Tr2i{Tn:£o'TF(^,y)} (26) 
\Pf{x)\' = Tr{TUTro'T^iHx,x)} (27) 

As in going from Eq. ^Jto Eg . 1131 the expectation values of the correlators can now 
be written in a form suitable for the implementation of the multi-level algorithm. In 
the expectation value of the Polyakov loop correlators one may take the average value 
{.)s inside a slab of each two-line operator appearing in Eq. ^J— EHHSZI To exploit 
the factorization into blocks along the direction 1, we evaluate the slab-averages as in 
Eq. El by taking the block average values of the line operators on the right-hand side of 
Eq.IHand Eq.^Hl 

It should be clear that the Polyakov loop correlator can be obtained in a factorized 
form for any desired representation. Also, the formulae presented hold for any number 
of dimensions d > 2. We now turn to some practical considerations concerning the 
implementation of the algorithm. 
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1. Produce an independent configuration by doing a sufficient number of hybrid over- 
relaxation sweeps through the whole lattice; 

2. for all slabs, in increasing time order: 



3. repeat the following procedure rus times: 

4. update the slab and then, for all blocks: 

5. do nib measurements of the needed line operators 
inside the block separated by an update within the block 

6. compute the two-line operators and add them to their slab averages; 

7. multiply the product of previous two-link operators by the new one; 



8. compute the traces. 



Table 1: Outline of the three- level algorithm. 







T/a 


L/a 


no 


ni 


n_L 


'^meas 


ro/a 


A20 


38 


20 


22 


4 


2 


11 


2536 


4.03879(39) 


A30 


38 


30 


22 


6 


2 


11 


2275 


4.03943(35) 


B 


44 


25 


26 


5 


2 


13 


1731 


4.82879(58) 


C 


54 


30 


30 


5 


3 


10 


1647 


6.1323(11) 



Table 2: The d = 3 SU(5) simulation details. 



2.2.3 Algorithm implementation and efficiency 

The three-level algorithm described above that we applied to d = 3 SU(5) simulations is 
in essence the Liischer-Weisz algorithm, where however the slab-averages of the two-line 
operators are computed in a factorized way. Due to the factorization along the direc- 
tion 1, we measure the Polyakov loop correlator in that direction only. The additional 
factorization breaks translational invariance and the symmetry between directions 1 
and 2 and thus leads to a loss of statistics; on the other hand the short-distance fluc- 
tuations of the one-line operators can be averaged out separately. Test studies jHE] in 
the d = 3 SU(2) theory show that the three-level algorithm improves on the original 
Liischer-Weisz algorithm for r > 0.9fm, a = 0.08fm, the latter being superior at shorter 
distances. This test concerns the static force in the fundamental representation, but in 
the higher representations, and also at smaller lattice spacing, we expect the additional 
factorization to pay off already at shorter distances. Since the correlation is obtained at 
fewer points, but each of them is evaluated more accurately, the memory requirement is 
lowered. 

The algorithm proceeds as sketched in Tab. ^ (the scope of loops are given by the 
indentation of the lines). The simulation parameters are given in Tab. El The number 
no of slabs and ni of blocks per slab can be found there. Because of memory limitations, 
we only computed the line-operators every second or third point in the 2 direction. The 
number n± in Tab.Elis the number of transverse coordinates at which the line-operators 
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were computed. We measured the correlation from distance 2a to 10a in all simulations; 
hence = 9. And in all simulations = 20 and = 60 (see Tab. ^ for the meaning 
of these parameters). We did not use the multihit technique [SZI, because for the higher 
representations it increases the number of two-line operator multiplications, which are 
expensive at large N. Finally, the total number of update sweeps between measurements 
is 60. They are grouped in compound sweeps consisting of one heat-bath followed by 
three over-relaxation sweeps [311 IS1 1^ • This proved sufficient to decorrelate almost 
entirely successive measurements of the observables and still represents a negligible over- 
head with respect to the measurements. To update the slabs we used three compound 
sweeps, and only one between line-operator measurements inside the blocks. 

The program goes from slab to slab in a time-ordered way. It is in order to save 
memory that each measurement of the two-line operator slab-averages is followed by 
their multiplication with the product of the preceding slab- averages. In this way, the 
program only needs enough memory for two fields of the size 

ni X n± X Ur X size(T). 

The two-line operator (|18|) contains A^^ complex numbers. Hence we reach the require- 
ment ~ 4 • 10^ 32-bit floating numbers for simulation C. 

For a given separation r, we correlate the line operators as far as possible from the 
boundaries. That is, for the correlator at a separation of 2na, we correlate the line 
operators computed at xi = Ai — na with those computed at xi = Ai + na (and all 
equivalent points by translation of distance Ai in the direction 1). For a separation of 
(2n + l)a, we average the correlation obtained at (Ai — na, Ai -|- (n + l)a) with that 
obtained at (Ai — (n + l)a, Ai + na). 

We introduce the notation 

rR(r) = (P^(x)PR(x + ri)) (28) 

for the on-axis Polyakov loop correlator. Tab.Olgives these correlators in the fundamen- 
tal, as well as in the AA-ality 2 representations 2a and 2s in simulation B. It turns out 
that these measurements in the different representations are quite strongly correlated 
at short distances (and essentially uncorrelated beyond r = tq), so that it will be ad- 
vantageous to take the ratio of static forces measured in the same simulation. With the 
tuning of the algorithm used, we are able to maintain a signal to noise ratio below 1% 
beyond Ifm in the fundamental representation and beyond 0.7fm in the 2a representa- 
tion. The performance (and memory requirements) of the algorithm behave favourably 
if the time extent is increased, as a comparison of the data from simulations A20 and 
A30 will reveal. 

To give an idea of the computing effort involved, one full measurement in simulation 
B took 1.7 hour on a single 2.4GIIz Opteron with 4GB of memory, while the 60 full- 
lattice sweeps between each of these full measurements took 61 seconds. 

The statistical errors are estimated with the Gamma method |42j . Simulation C 
was run with 8 replica (= independent simulations) and the others with 4 replica; we 
checked in all cases that the results are consistent across replica. We also found that 
using the jacknife method gave error estimates in agreement with the Gamma method, 
which is hardly surprising since the successive measurements were very decorrelated in 
these runs. 
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3 The static force in the fundamental representation 



In this section we present our data on the static force. It is found that the quantity TqC" 
has very small discretization errors. By comparison with the SU(3) data of Liischer and 
Weisz jin] and the SU(2) data of Majumdar |3T| (see also p2J), we find that it is also 
very weakly dependent on the number of colors A^. We then compare the N = 2,3,5 
results for the effective central charge at fixed value of rg. 

3.1 The SU(5) static force 

The static force is defined at finite lattice spacing by 

Fn{r) = ^[VK{r + a)-VR{r)], Vn{r) = -^logTn{r). (29) 

We reuse the tree-level improved argument f = r + a/2 + 0(a^) of the function given 
in ^U]. The fundamental static force obtained in simulation B can be found in Tab.|lJ 
The reference scale ro/a given in Tab. [21 is obtained in the standard manner ^2]; by 
linearly interpolating (r/a)'^ as a function of r'^F{r) to the abscissa 1.65. 

To a first approximation, the force is constant beyond 0.75fm, signalling linear con- 
finement. To investigate this in more detail, we find it useful to define the quantity 

^W = i^W-^- (30) 

If the Liischer term does provide the leading asymptotic correction to the linear poten- 
tial [S] , then this quantity deserves the name of effective string tension and the notation 
is appropriate. Note that by definition, 

r2CT(ro) = 1.65 - ^ = 1.519100. . . (31) 
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at all lattice spacings. In particular, discretization errors on rQa{r) are automatically 
reduced around separation r = tq. 

The function rQcr(r) is plotted in Fig. I^lfor all data sets. The practical advantage of 
considering this function is that it is very flat beyond 0.4fm (note the scale of the vertical 
axis); in this way, the data can be visualized in more detail. We firstly remark that the 
data A20 and ^30 agree within error bars at all separations; the common abscissa of 
these data points have been split symmetrically for better visibility. This indicates that 
the time extent T ~ 5ro is sufficient to 'filter out' the ground state of the transfer matrix 
in the presence of the two Polyakov loops ^01 separated by r < 2ro; the statement holds 
of course at the precision to which we have been able to compute the static force. 

Secondly, the scaling violations of rQa{r) are at the few permille level (or less) around 
r = IAvq. To carry out a continuum extrapolation, we need to interpolate rQF(xro) to a 
few fixed values of x. This is done by a linear interpolation in 1/x^ (in 10^ a three-point 
interpolation in 1/x^ was used). The result is given for each simulation in Tab. IHl It 
turns out that extrapolating these quantities to the continuum with the expected O(a^) 
corrections jl3] yields poor values (see Fig. El). Therefore we extrapolate the data 
from the two smaller lattice spacings linearly in (a/rg)^, and use a quadratic fit to all 
three lattice spacings to estimate the uncertainty associated with the linear fit; we refer 
the reader to the appendix for more details. At least one additional simulation at a 
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Figure 2: The SU(5) effective string tension in units of Tq ^. By interpolation all data 
go through the cross (by definition of ro). 

lattice spacing of a ~ O.lro is required to take the continuum limit in a way that does 
not blow up the final uncertainty. 

Finally, we did not investigate finite- volume effects, but since 2.45fm < L < 2.72fm 
in these simulations, we do not expect them to affect our conclusions. 

3.2 The quantity arg for all 

In comparing SU(A^) gauge theories, perhaps the most natural way to set the common 
energy scale is by equating their (fundamental) string tensions a. This quantity is 
however more difficult to compute with high precision than tq, because it involves an 
extrapolation to infinite string length |15j. For that reason, we are now going to show 
that the dimensionless quantity TqU is A'^-independent to a precision better than 1%. 
As a consequence, comparing SU(A^) gauge theories at fixed value of rg is essentially 
equivalent to fixing the value of a. 

We can extract the effective string tension at Ifm straightforwardly. A direct com- 
parison of our data with the SU(3) data of Liischer and Weisz at the same lattice spacing 
a = 0.207ro gives 

rgcT(2ro) = 1.5124(65), SU(5), /3 = 44, (32) 
r^a{2ro) = 1.51776(24), SU(3), (3 = 15. (33) 

Interestingly, the string tension has also been computed in the closed string channel by 
Teper |4j at these values of (3. Using the values of ro quoted above, we have 

SU(5) : r2 adosed = 1.521(8) (34) 
SU(3) : r2 aeioscd = 1.520(14). (35) 



10 



0.16 



Effective central charge 



-1 1 1 1 1 r- 



5= 

o 



0.14 
0.12 
0.1 
0.08 
0.06 



p=54 
p=44 
T/a=30 p=38 
T/a=20 B^38 



_i L I I L_ 



_l L_ 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 
r/rn 



Figure 3: The SU(5) effective central charge at several lattice spacings. The perturbative 
prediction is shown at short distance, and the horizontal line marks the value 

The string tension quoted by Teper is also an effective string tension (as defined by 
m{L) = a{L)L — but at a closed string length of 3.3ro, where higher order string 
corrections are presumably much reduced. In the SU(3) case an additional calculation of 
his at string length of S.Oro yielded the same string tension (and accuracy). We conclude 
that the closed string tension extracted from torelon masses agrees with the open string 
tension extracted from the static force at Ifm at a precision of 1%. 

Majumdar [SJI reports almost identical values of r^a in the SU(2) theory at lattice 
spacings a < O.lGrQ. In particular, a comparison between SU(2) and SU(3) closer to the 
continuum yields 

rla{2ro) = 1.523(2) SU(2), /3 = 7.5, ro = 6.29a (36) 

rgf7(2ro) = 1.51999(34) SU(3), /? = 20, ro = 6.71a. (37) 

Cutoff effects on TqCT are thus seen to lie well below the percent level. Assuming that 
no strong A^-dependence appears beyond N = 5 (the contrary would be very surprising 
in view of the results obtained in Q), we reach the following remarkable conclusion: 

r^a = 1.52 • (1 ± 1%) ViV > 2, a < 0.21ro. (38) 

In particular, 

a = CT(ro) • (1 ± 1%) ViV > 2 (39) 

holds in the continuum limit (see also Fig. ^ . 

3.3 The effective central charge 

The effective central charge is defined by 

c(f) = -if3[1/(r - a) - 2V{r) + V{r + a)]/a^, (40) 
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where f = r + 0(a^) is defined in TO^. Fig.|31shows the computed effective central charge 
c(r) in the range of distances 0.25 to 0.75fm. A precocious convergence to values within 
5% of ^ is observed. 

By eye, the data at /3 = 44 and /? = 54 practically fall on top of each other up to 
0.6fm, indicating that discretization errors are reasonably small. Beyond that point, 
the /? = 54 data's accuracy no longer allows for a useful comparison. The data points 
obtained at the coarsest lattice spacing systematically lie slightly above the other points. 
This effect was already seen in SU(2) and SU(3) calculations. Also, the data sets ^420 
and A^Q are well compatible, indicating that the time extent was chosen long enough 
for the present statistical accuracy. 

The expected behaviour at short distances, as given by two- loop perturbation the- 
ory [22], is shown on the figure, as well as the universal value ^ expected at long dis- 
tances. The continuum value of g'^rQ, which determines the slope of the short-distance 
prediction, was estimated by extrapolating the value of g'^r^ of simulations B and C 
linearly in g'^a. Using the continuum value of \fo j g^ of 4 and Eq. LSSI vields a slope 
smaller by 2%; the difference would hardly be noticeable on Fig. |21 

In view of taking the continuum limit, we interpolate cixro) to a handful of values 
of X. We use a three-point polynomial interpolation in 1/x (following ^3); the result 
is given in Tab. d The continuum limit is described in the appendix, where also an 
alternative definition of the effective central charge is proposed and shown to lead to 
the same continuum results (Fig. |H1). 

3.4 The effective central charge for all N 

It is now interesting to compare the effective central charge curves obtained for = 2, 3 
and 5. If we attempted to compare the data in the continuum limit, the continuum 
extrapolation would lead to the dominant source of uncertainty, because only three 
lattice spacings are available at each A^. Also, the data at the coarsest lattice spacing in 
the SU(3) and SU(5) data sets are not (always) consistent with O(a^) corrections. For 
that reason we prefer to compare the data sets at the smallest common lattice spacing. 
This corresponds to /3 = 7.5, 20 and 54 for A^ = 2, 3 and 5 respectively. The values of 
ro/a, 6.29, 6.71 and 6.13 respectively, indeed are matched at the 10% level. 

The comparison is illustrated on Fig.|^ At short distances, the curves are expected to 
differ by 0(1/A^^) from perturbation theory and the scaling of g'^Nr^. This is indicated 
by the straight solid lines which stop short at 0.4ro around where perturbation theory 
ceases to be accurate. Remarkably, all three curves seem to approach the value of ^ 
at larger distances. This fact constitutes strong evidence for the universality of the 
string correction in d = 3 SU(A^) gauge theories. On the other hand, the approach 
to the asymptotic value is surprisingly precocious. The numerical closeness of c(r) 
to ^ is probably deceptive, in so far that it suggests that any higher order string 
corrections are already very small below Ifm. The leading correction term is 0(l/cTr^), 
and its coefficient is universal and positive in the effective theory describing the massless 
transverse fluctuations of the flux-tube jlOj . The Nambu-Goto string in a non-critical 
number of dimensions corresponds to a special case in this class of theories; it provides a 
consistent analytic prediction for the r x T Polyakov loop correlator to any finite order 
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Figure 4: The effective central cfiarge for several gauge groups. The SU(3) data is taken 
from |10j, the SU(2) data from ^j. The perturbative prediction is shown at short 
distance (the slope increases slightly with A^) and the Nambu-Goto prediction at long 
distances, while the horizontal line marks the value 

in 1/r or 1/T |lUj . It predicts an effective central charge [l^ of the form 

^'^^'M-^ (l-T^T.)"''. (11) 

shown for illustration in the upper right corner of Fig. \^ It suggests that the effective 
string theory becomes accurate at significantly larger distances. To summarize, since 
c(r) < ^ at short distances; since the 1/r series for c(r) at large distances is almost 
certainly asymptotic; and since the 0(l/cjr^) term is positive, this term improves on 
the leading order prediction at the earliest when the effective central charge computed 
in the gauge theory rises above the value 

The preceding arguments thus motivate the definition of the distance rg where the 
effective central charge crosses ^: 

cirs) ^ f,. (42) 

This distance is a natural separation scale between the perturbative regime and the string 
regime. The ratio Tg/ro is expected to have a continuum limit with O(a^) discretization 
errors, and at finite lattice spacing an interpolation formula must be prescribed to define 
it. In general we propose to use a three-point polynomial interpolation in r using the 
three nearest data points, because c(r) has quite some curvature around rg. 

Fig. |l]shows quite strikingly that as the number of colors is increased, the curve c(r) 
interpolates over an even shorter distance range between the perturbative behaviour and 
the value We therefore study the A'"-dependence of rg/rQ. To do so we shall presently 
use a simple linear interpolation between the two nearest points to the crossing, because 
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d=3 SU(N) LGT (a/ro~0.16) 




Figure 5: The dependence of rg/ro on the number of colors N. The SU(3) data is taken 
from the SU(2) data from 



the SU(2) and SU(5) data rapidly become less accurate beyond the crossing. Luckily in 
all three cases there is a point lying almost exactly at which makes the details of the 
interpolation a negligible source of uncertainty. In order to have the same discretization 
scheme for all data sets, we replaced r by f in the SU(2) values of the effective central 
charge We thus find: 



1.79(5) N = 2 
1.32(5) N = 3 
0.96(2) N = 5 



(43) 



The smaller uncertainty at SU(5) reflects the fact that c(r) cuts the horizontal ^ line at 
a shorter distance and with a larger slope than for SU(2) and SU(3). At = 5, /? = 44 
we obtain Vs/vq = 0.977(11), in agreement with (|43j) . The trend that rg/ro strongly 
decreases with the number of colors is quite striking and constitutes a rare example of 
a physical quantity in SU(A) gauge theories which has a strong A^-dependence. This 
dependence is shown on Fig. The data is well compatible with a linear function 
in 1/A: 

^ = 0.402(45) + 2.80(20)-^ {x^/d.ol ~ 0.01). (44) 

Could it be that rg becomes as small as 0.2fm in the large-A limit? An SU(8) calculation 
would provide a useful clue to answer the question. In the mean time, we can compare 
the value to two other distance scales which are natural in this context. The first comes 
from the Arvis formula, Eq. 1411 The latter becomes singular at a distance given by 



— — = 0.415(4). 
120- ro ^ ^ 



(45) 
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3D SU(5) gauge theory 



0.005 



p=54 
(3=44 
p=38 T=30 
B=38 T=20 




Figure 6: Relative deviation of the force ratios F^/Fp from Casimir scaling in the SU(5) 
gauge theory. The points above and below zero correspond to the representations R=2a 
and R=2s respectively. 

where we have used Eg. 1381 It is intriguing that the naively obtained value rs{N = oo) 
is compatible with r^. Secondly, the dotted line on Fig. |3] is the continuation of the 
curve corresponding to the perturbative two-loop result j29( in the large- limit (using 
limAT^oo V^/g'^N = 0.1975(10) from |lj and Eq.EHI)- It cuts ^ at distance rpt given by 

^ = 0.527(4). (46) 

The rise of c(r) up to and beyond which becomes more rapid as N is increased, 
could be a sign that the asymptotic series in [ar'^)~^ provided by the effective string 
theory converges better at a given value of ^/o^r if the number of colors is increased. 
This is not unconceivable, because certain effects such as the radiation of glueballs by a 
highly excited string, which exist in the gauge theory but are neglected in the effective 
string theory, are l/N'^ suppressed. 

4 What is the size of Casimir scahng violations? 

We now address this question in the case of the two smallest irreducible representations of 
SU(5) of AA-ality 2. These are obtained by taking the direct product of two fundamental 
charges and forming the (anti) symmetric linear combinations, as in Eq. [3 and IHl We 
discuss the ratio of these forces to the force in the fundamental representation. Note 
that, unlike the ratio of static potentials, this quantity has a continuum limit. 
The Casimir scaling prediction is, in a range of distances r^iri < r < r^ax, 

F2.{r) ^ 2{N - 2) 

Fi(r) Ci A^-1 



3 
2' 



(47) 
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F2s{r) cj C2S ^ 2{N + 2) ^7 

Fi(r) Ci iV + 1 ;v=5 3- ^ ' 

Perturbation theory predicts j^H] 

FR{r) = ^ [l + l^g^Nr + ©((^^iVr)^)] . (49) 

In particular, Casimir scaling of static forces holds at least to order g^. 

The result for the ratio of static forces obtained in simulation B is given in Tab. 0J 
Fig. ini shows the Monte-Carlo data on the relative deviations of i^2a/-^i from Casimir 
scaling in the SU(5) lattice gauge theory. The window in r is 0.2 — 0.7fm. Notice the 
scale of the vertical axis: the (2a) data points indicate a positive relative deviation which 
grows from about 0.001 to 0.01 in that range of distances. The (2s) data on the other 
hand shows a stronger, negative deviation. This is expected, since the static charges in 
this representation can be screened to (2a). The force between charges in the adjoint 
representation, which can be screened completely, also exhibits this behaviour [2^]. A 
positive deviation is thus rather special. It will remain positive at larger distances if the 
k = 2 string tension is larger than C2a/ Ci ■ or if it is given exactly by that expression 
and its central charge by 

Beyond ro, we are perhaps seeing the effects of finite time extent: there is a 1.9a 
discrepancy between the third (2a) data points of simulations A20 and A^q on Fig. El If 
present, this effect comes from F2a., because the quantity rQFi{r) is perfectly consistent 
between A20 and ^30 (see Tab. El- 

It is clear from Fig. IHlthat the (absolute) discretization errors are bounded by 0.001. 
The force ratio -F2a/-^i at distance xro has been interpolated to a few values of x in 
Tab. [TUl The two-parameter function f{x) = px'^/{xP' + q^) was used to interpolate 
the quantity plotted on Fig. El but the result never actually differs by more than one 
standard deviation from a simple linear interpolation in x. Tab. 1101 also gives the result 
of a continuum extrapolation performed along the same lines as for rgFi(r). 

Since there is a stable string in the AA-ality 2 sector, it makes sense to consider the 
effective central charge in the 2a representation; Fig. [7| shows our data on C2a('^)- At 
short distances, it is roughly a factor 3/2 larger than in the fundamental representation, 
as predicted by perturbation theory. The curve seems to flatten off, as it does in the 
fundamental representation, but the data does not extend far enough to show whether 
it will bend down towards ^, as we expect it to jlSl. It was observed (30j in closed 
k = 2 SU(8) strings that C2a < ^ci- 

4.1 Mixed representation correlators 

It is interesting to look at the cross correlation {P2a.{x)* P2s{x + fl))- There is no sym- 
metry which makes this correlator vanish. However we find that it is consistent with 
zero. For instance, the 'overlap' 



(P2a(x)*P2s(^ + ri)) 

v/r2a(r)r2s(r) 
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Figure 7: The SU(5) effective central charge in the 2a representation. The perturbative 
prediction is shown at short distance and the Nambu-Goto prediction at long distance. 
The horizontal line marks the value tSt. 



is consistent with zero in simulation B with a statistical uncertainty of 5 • 10 ^,3-10 ^ 
and 4 • 10~^ at distances r = 2a, 3a and 4a respectively. Since 

{P2.(xr P2s{y)) = i {{{Tt {U*p})\x) {TT{Up})\y)) - (Tr{C/2}*(x) Tv{Uj,}{y))) , 

the statement that the cross correlation vanishes is equivalent to the property that the 
two correlators on the RHS of this equation are equal. 

The interpretation that we propose is based on the observation that the direct pro- 
duct of the representations (2a)* and 2s does not contain the singlet representation 
(though it does contain the adjoint representation). Therefore if we imagine inserting 
the static particles adiabatically into the system, they initially have an overall color 
charge. The only way the free energy can be made finite is if virtual adjoint charges 
screen the system. Since the external color charge object is extended, it presumably 
costs a large amount of energy to screen it. In addition, statistically speaking the color 
of the virtual adjoint charge must match that of the static system, which happens with 
a probability of order l/N"^. Thus one expects the overlap to be l/N"^ color-phase- 
space suppressed j2B]- All this results in a very small cross-correlation of the Polyakov 
loops. This large- A'^ argument was already made in [30] to explain the fact that in 
closed /c-string computations |28[ I3()j . the wavefunction of the lightest state is very well 
approximated by a (fuzzy) Polyakov loop in the totally antisymmetric representation of 
AA-ality k. 
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5 Conclusion 



We computed the static force in the d = 3 SU(5) gauge theory employing an efficient 
three-level algorithm. Linear confinement is observed, and the effective string tension 
defined in Eq. IHillis essentially constant beyond 0.5fm. The string tension extracted at 
Ifm agrees with the closed string tension extracted from torelon spectroscopy jl]. 

By comparison with existing SU(2) [221331 and SU(3) ^01 data, the quantity ar^ is 
found to be independent of the number of colors at the 1% level; it also has very small 
discretization errors. Thus comparing SU(A^) gauge theories at common string tension 
or common Sommer reference scale is equivalent at that level of precision. 

The effective central charge c(r) was obtained in the range 0.25 — 0.75fm. It converges 
to within 5% of the expected asymptotic value of confirming the multiplicity and 
the bosonic nature of the fiux-tube's massless degrees of freedom. A comparison with 
the SU(2) and SU(3) data reveals that the distance rg where c(r) crosses the value ^ 
decreases steadily by almost a factor two from SU(2) to SU(5). Since (bosonic) effective 
string theory predicts that the asymptotic value is approached from above |31, it is 
tempting to speculate that the asymptotic expansion's accuracy is higher at a fixed 
value of ^/o^r for larger N. 

We also studied the static force in the symmetrized (2s) and antisymmetrized (2a) 
direct product representations of two 'quarks'. For > 4, it is known [2H] that such 
static charges of AA-ality 2 are linearly confined with a different string tension from 
the fundamental one. And indeed we find that the ratio F2a,/Fi is constant to a first 
approximation, confirming the linear confinement property. Furthermore, it is very close 
to the Casimir scaling prediction (3/2 in this case). A more detailed study reveals that 
deviations are present at the 0.1% to 1% level in the range 0.2 — 0.7fm. They are positive, 
unlike the case of the adjoint static force where screening of the adjoint fiux must 
set in at some distance r. The 2s representation on the other hand exhibits a negative 
deviation, as one expects if there is a single stable string per AA-ality and screening 
occurs. 

Since the representation 2a has AA-ality 2 and has the smallest quadratic Casimir of 
all irreducible representations in that sector, one might expect to find again a precocious 
onset of the effective central charge C2a{r) for the k = 2 string. This is not the case: at 
r = ro, where c(r) lies within a few percent of C2a{r) is almost as large as | • 
Data at further distances is needed to see whether it decreases towards 
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Appendix: discretization errors and the continuum limit 

The quantity rQF(r) has been extrapolated hnearly in {a/rof to the continuum, using 
the data from simulations B and C. This is our best estimate in the continuum limit 
and it can be found in Tab. IHl The first number in brackets is the statistical error on 
the result (there is no x^)- To estimate the systematic uncertainty on this procedure, 
we also do a quadratic fit in (a/ro)^ to simulations ^30, B and C (again, there is no 
X^). The second number given in brackets is the continuum value obtained by quadratic 
extrapolation minus the quoted continuum result. Thus a conservative error estimate 
is the maximum of the absolute value of the two numbers in brackets. The continuum 
limit is illustrated for rQa(r), which is equivalent to rgF(r) for this purpose, on Fig. |H1 

We proceeded in the same way for F2a{f') / Fi{r) (Tab. [T?H) . For the effective central 
charge, we also give the continuum c(r) obtained in this way in Tab. d Recall that the 
definition of Ref. was used to define c(r) at finite lattice spacing. We now show that 
a different definition gives the same result in the continuum limit. 

We propose to determine the parameters a?a{r), afl{r) and c(r) by 

c(r) = + ^v{r-a)-2V{r) + V{r + a)) (50) 

a'^a{r) = ^{r-a)V{r-a)-rV{r) + l{r + a)Vir + a) (51) 

a~,{r) = "^(^ + V{r - a) + V{r) + § ) y(, + (52) 

a a a 

This new definition is such that if the on-axis lattice Polyakov loop correlator was given 
by 

g-Vb{r)r^ Fo(r-) =/i + (7r--, (53) 

r 

then the effective parameters match those appearing in Vo(r) exactly. Note that the 
functions c and c differ only by 0{a?) terms, and that c(r) is simply a factor (1 — {a/r)"^) 
times the 'naive' definition of the effective central charge used in j31j . And indeed, 
comparing Tab. [H] to Tab.[7J we find that c(r) is compatible with c(r) in the continuum 
limit, as it should be; see Fig.|Hlfor an illustration. This consistency-check is also a test 
of our error analysis. 
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Effective string tension: continuum extrapolation 
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Figure 8: Continuum extrapolation of the effective string tension r^a^r) = rQF{r) — 
^(ro/r)^ and two definitions of the central charge c(r) and c(r) in d = 3 SU(5) gauge 
theory. The result, with its final error bar, is shown at vanishing abscissa. 
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Table 3: The raw data of simulation B. 
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F2a{r)/Fi{r) 


F2s{r)/Fi{r) 




2.379 


0.0857165(97) 


1.501928(46) 


2.32216(20) 




3.407 


0.076183(12) 


1.504059(99) 


2.3092(28) 




4.432 


0.071825(15) 


1.50613(20) 


2.424(95) 




5.448 


0.069546(18) 


1.50793(71) 






6.458 


0.068237(23) 


1.5144(33) 






7A6A 


0.067463(39) 


1.507(33) 






8.469 


0.066812(67) 








9.473 


0.06634(25) 







Table 4: The static forces in simulation B. 



r/a 


Cl(f) 


C2a{f) 


C2s{f) 


2.808 


0.105535(54) 


0.15671(10) 


0.2560(23) 


3.838 


0.12320(31) 


0.18110(61) 




4.875 


0.13202(48) 


0.1916(28) 




5.902 


0.1346(21) 


0.158(24) 




6.920 


0.1283(45) 






7.932 


0.162(16) 







Table 5: The effective central charges in simulation B. 
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B 
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1.97861(46) 


— 


0.65 


1.81564(22) 


1.81552(18) 


1.81409(28) 


1.81380(35) 


1.8133(10)(20) 


0.80 


1.72229(22) 


1.72211(18) 


1.72119(15) 
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1.7209(8)(14) 


0.95 


1.664406(60) 


1.664354(48) 


1.664300(52) 


1.66417(12) 


1.6640(3)(-2) 


1.00 


1.65 


1.65 


1.65 


1.65 


1.65 


1.05 


1.637603(52) 


1.637647(41) 


1.637694(45) 


1.63780(10) 


1.6380(3)(1) 


1.20 


1.60865(24) 


1.60849(16) 


1.60934(19) 


1.61001(44) 


1.6111(12)(-2) 


1 Q r; 


i.0OO00^4Dj 


i.oooio^oU j 


1 ^Q.Cill ( A'7\ 

l.ooy ( (\^() 






1.50 


1.57272(81) 


1.57386(65) 


1.57638(76) 


1.5686(32) 


1.556(9)(-20) 


1.65 


1.5615(18) 


1.5635(19) 


1.5647(10) 






1.80 


1.5536(47) 


1.5558(55) 


1.5551(23) 






2.00 


1.562(13) 


1.548(32) 


1.5451(65) 







Table 6: The quantity r^Fir) interpolated, and its continuum limit (details in the App.). 



r/rQ 


^20 


^30 


B 


c 


cont. 


0.458 








0.090714(40) 




0.65 


0.11028(74) 


0.11141(56) 


0.11284(27) 


0.11163(19) 


0.1096(7)(-50) 


0.80 


0.12573(64) 


0.12484(50) 


0.12348(31) 


0.12202(28) 


0.1196(9)(-5) 


0.95 


0.13282(56) 


0.13235(44) 


0.13006(36) 


0.1304(10) 


0.131(3)(5) 


1.00 


0.13424(57) 


0.13412(42) 


0.13172(46) 


0.1314(11) 


0.131(3)(4) 


1.05 


0.13534(71) 


0.13564(48) 


0.13279(44) 


0.1315(12) 


0.129(3)(2) 


1.20 


o.i;]74(;i.')) 


0.1390(10j 


0.1344(19) 


0.l;]()2(31j 


0.139(8)(;i2) 



Table 7: The effective central charge c(r) interpolated, and its continuum limit. 



r/ro 


-420 


-430 


B 


C 


cont. 


0.458 








0.098331(44) 




0.65 


0.1105(19) 


0.1135(15) 


0.11787(15) 


0.11611(21) 


0.1132(6)(-114) 


0.80 


0.13196(48) 


0.13130(37) 


0.12939(35) 


0.12533(27) 


0.1187(9)(-48) 


0.95 


0.13995(64) 


0.13921(51) 


0.13491(33) 


0.13260(92) 


0.129(2)(3) 


1.00 


0.14117(58) 


0.14074(46) 


0.13606(42) 


0.1340(11) 


0.131(3)(4) 


1.05 


0.14194(62) 


0.14190(45) 


0.13690(45) 


0.1339(12) 


0.1289(32)(7) 


1.20 


0.1425(14) 


0.14393(92) 


0.1375(16) 


0.1367(27) 


0.135(7)(6) 



Table 8: The effective central charge c(r) interpolated, and its continuum limit. 
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T It(\ 


B 


c 


0.458 




0.135447(60) 


0.65 


0.16725(88) 


0.16486(32) 


0.80 


0.18146(60) 


0.17752(69) 


0.95 


0.1894(20) 


0.1849(30) 


1.00 


0.1912(27) 


0.1860(54) 


1.05 


0.1868(38) 


0.186(12) 


1.20 


0.162(21) 


0.135(51) 



Table 9: The quantity C2a(r) interpolated. 



r/ro 


A20 


^30 


B 


C 


cont. 


0.50 






1.324(31) 


1.130(39) 




0.65 


2.434(65) 


2.567(50) 


2.279(50) 


2.087(48) 


1.77(15)(13) 


0.80 


3.59(15) 


3.693(97) 


3.317(66) 


3.32(12) 


3.3(3)(7) 


0.95 


5.09(36) 


4.01(61) 


4.28(15) 


4.80(21) 


5.7(6)(6) 


1.00 


5.66(67) 


3.94(81) 


4.57(21) 


5.35(29) 


6.6(8)(4) 


1.05 


6.25(84) 


3.87(97) 


4.86(30) 


5.92(43) 




1.20 




8.5(3.2) 


6.57(59) 


7.0(2.2) 





Table 10: The quantity 10^ • {^F2a{r) / Fi{r) — 1) interpolated, and its continuum limit. 



r/ro B 


C 


0.50 -4.943(88) 


-4.55(10) 


0.65 -8.68(69) 


-8.73(57) 



Table 11: The quantity 10^ • (f F2^(r)/Fi(r) - 1) interpolated. 
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